Influence of polarization of laser beam on emission intensity of femtosecond laser-induced breakdown spectroscopy
Yang Lan1, 2, Liu Miao1, Liu Yi-Tong1, Li Qing-Xue1, Li Su-Yu1, Jiang Yuan-Fei1, Chen An-Min1, †, Jin Ming-Xing1
Institute of Atomic and Molecular Physics, Jilin University, Changchun 130012, China
College of Physics, Jilin University, Changchun 130012, China

 

† Corresponding author. E-mail: amchen@jlu.edu.cn

Project supported by the National Natural Science Foundation of China (Grant Nos. 11674128, 11674124, and 11974138), the Jilin Provincial Scientific and Technological Development Program, China (Grant No. 20170101063JC), and the Scientific and Technological Research Project of the Education Department of Jilin Province, China (Grant No. JJKH20200937KJ).

Abstract

Laser-induced breakdown spectroscopy (LIBS) is an important technique which is widely used to analyze element composition. In order to improve the sensitivity of LIBS, much effort has been made to enhance the spectral intensity of LIBS by proposing a number of methods. In addition, we find that laser polarization has great influence on the emission intensity of femtosecond LIBS. By comparing the emission intensity of femtosecond LIBS in the circular polarization with that in the linear polarization, the spectral intensity in the case of circular polarization is stronger than that in the case of linear polarization. Moreover, this phenomenon is more obvious as laser energy increases. The polarization plays an important role in LIBS signal intensity. Based on the observation, the enhanced mechanism of the laser polarization for the spectral intensity is discussed in this paper, which will be helpful in spectral analysis and component analysis.

1. Introduction

Laser-induced breakdown spectroscopy (LIBS) is an effective technique to analyze the component of material without any preparation of sample.[15] The principle of this method is as follows. A high-energy laser beam is focused onto the target surface to produce plasma, and then plasma creates the emission spectrum: this emission spectrum is used to analyze the element composition.[69] In recent years, the LIBS has achieved dramatic development and received much attention due to its practical applications such as in environmental monitoring, biological agents, and aerosol testing. More importantly, the high efficacy, non-intrusiveness, and real-time analysis make it extensively employed for detecting the sample composition. Although the LIBS has been successful in numerous applications, its sensitivity, accuracy, and precision are still poor compared with those of other analytical techniques.[10] Through improving the spectral intensity, the sensitivity can be greatly enhanced since the analysis process is easier and clear. In the past decades, many people have proposed some techniques to improve the spectral emission intensity, such as double-pulse excitation,[1015] electromagnetic field,[1620] spatially confined LIBS,[2125] nanoparticle-enhanced LIBS,[26,27] femtosecond pulse excitation,[28,29] spark-assisted LIBS,[3033] and increasing sample temperatures.[3438] In addition, changing the polarization state of laser pulse may also be an effective and feasible method, which can be easily achieved by using a quarter-wave plate.

The phenomenon that the spatial distribution of the electric vector vibration of light wave loses symmetry in the direction of light propagation is called the polarization of light, and only the shear wave can produce polarization. Any light whose vibration loses this symmetry is collectively called polarized light. The polarization states can be divided into three types: the circular one, the linear one, and the elliptical one. The polarization characteristics of laser pulses can be changed by rotating 1/4 wave plates: circularly polarized laser is produced when the rotation angle of 1/4 wave plates is 45° + n × 90°, and linearly polarized laser is created when the rotation angle is n× 90°. Polarization state of laser pulse can affect the spectral intensity through influencing the electron density and plasma temperature, which are the most important characteristics of plasma.[39] We experimentally demonstrated that circularly polarized laser pulses can produce the spectra with higher emission intensity than the linear laser pulses. In addition, the width of laser pulse also acts as an essential feature in investigating the effect of laser polarization.[40]

With the fast development of laser technology, femtosecond pulse laser can output higher laser energy. The femtosecond pulse laser is accompanied with a higher intensity, hence there is a significant benefit in using ultra-short pulse laser: its output power is high enough to be applied to many fields.[4146] The ultra-short laser pulses have been used for examining and processing the materials in the past decades,[4749] and they can improve the ablation repetition rate and reduce the sample pollution. More importantly, the ablated hole size is small, enhancing the control of micro-fabrication process, when comparing with the scenario of the nanosecond pulse or longer pulse.[5054] Utilizing the advantages mentioned above, femtosecond laser as an excitation source is introduced for improving the sensitivity of the LIBS. Femtosecond laser with ultra-fast and ultra-intense characteristics can rapidly inject all energy into the sample based on its very short pulse width. As a result, the spectral emission produced by femtosecond LIBS has some advantages including a low continuous spectrum, allowing the spectral signals with many laser shots to accumulate.[28] However, the reduced interaction between excitation light source and produced plasma causes emission intensity to lower.

Referring to different laser polarization effects on LIBS, femtosecond laser is introduced into LIBS. Comparing with optical emission induced by the linearly polarized laser, and changing the polarization of laser pulse into elliptical and circular cases, the emission intensity of LIBS may be obviously improved. Recent investigations of electron momentum distribution after the tunnel has been ionized by femtosecond laser indicate that the partitioning of energy for circular polarization is different from that for linear polarization. For instance, Mitryukovskiy et al. studied that the optical emission from neutral and singly ionized N2 molecules created by femtosecond laser filamentation in atmospheric environment depends strongly on femtosecond laser polarization,[55] and their results illustrated that the femtosecond laser polarization states influence the plasma through affecting the kinetic energy of electrons. Lemos et al. observed that the expansion of the plasma plume created by circularly polarized laser is quicker, which indicates that the electron temperature is increased by more than 50%.[56] Wang et al. demonstrated that circularly polarized femtosecond laser can improve emission intensity of LIBS.[57] The polarization can play an important role in femtosecond LIBS signal intensity.[57] However, the physical reason of the influence of polarization state on the spectral intensity is less investigated and still is an open question. Therefore, we investigate the influence of the polarization on the spectral emission intensity in detail.

In this paper, we discuss the influence of the polarization of laser pulse on spectral intensity by comparing the breakdown spectra of Cu, produced by different polarized laser pulses. The spectral intensity values produced by femtosecond laser pulses under different polarization states are investigated, thereby trying to figure out the physical pattern of effects caused by polarization. Meanwhile, we also calculate the plasma temperature and electron density under different polarization states, which are both essential characteristics of plasma and can conduce to explaining the mechanism of this phenomenon.

2. Experimental setup

Figure 1 shows the schematic diagram of experimental setup. All the measurements were carried out in the atmospheric environment. The femtosecond laser system used is a one-box ultrafast Ti:sapphire amplifier (Coherent Libra) with a wavelength of 800 nm, width of 50 fs at a repetition rate of 10 Hz. And, the synchronization and delay generator (SDG) of the femtosecond laser system was used as a trigger source for the detection system. A half-wave plate and a Glan laser polarizer were combined to change the laser energy, and a quarter-wave plate was used to change the laser polarization. A focusing lens with 20-cm focal length focused the laser perpendicularly onto a target. The target (pure Cu, 99.9%) was attached to a three-dimensional translation stage. The produced plasma spectra were guided to a spectrometer (PI Action, 1200 lines per mm) through a fiber. The spectra were detected by using an ICCD (PI-MAX4) with 0-μs delay time and 10-μs gate width. The spectral resolution was 0.04 nm. Each spectrum was an average of 100 shots. In the present measurement, the laser polarization was adjusted by rotating the quarter-wave plate an angle. When its rotated angle is 45° + n × 90°, the laser polarization is circular; while the angle is n × 90°, the laser polarization is linear.

Fig. 1. Schematic diagram of experimental setup for Cu plasma produced by femtosecond laser for different polarizations (M: mirror, HWP: half-wave plate, Glan: Glan laser polarizer, Q: quarter-wave plate, L: lens, and ICCD: intensified CCD).
3. Results and discussion

Figure 2(a) shows the emission spectra of the plasma generated by femtosecond laser pulses with linear and circular polarizations, respectively. The laser energy is 1.9 mJ. From this figure, the spectral intensities of Cu (I) lines at the wavelengths of 510.55 nm, 515.32 nm, and 521.82 nm are stronger in the case of circular polarization than in the case of linear polarization. Figures 2(b) and 2(c) show the peak intensities of Cu (I) lines at 510.55 nm and 521.82 nm varying with laser energy. From these figures it follows that the spectral intensity for the circular polarization is higher than that for the linear polarization. When the pulse energy is 0.5 mJ, the spectral peaks tend to the same intensity for linear and circular polarization (The pulse energy at which spectra share the same intensity is called “energy threshold”). Continue to increase the laser energy, then the signals at 510.55 nm and 521.82 nm will become stronger for the circularly polarized laser. As the energy of the laser pulse increases, the phenomenon becomes more obvious.

Fig. 2. (a) Comparison between spectral intensities of linearly and circularly polarized femtosecond laser-induced Cu plasma with laser energy bing 1.9 mJ. Spectral peak intensities of Cu (I) lines at (b) 510.55 nm and (c) 521.82 nm varying with laser energy.

In femtosecond laser ablation of metal, the absorption mechanism is the mechanism for inverse bremsstrahlung absorption.[58] Femtosecond laser ablation of metal is a non-equilibrium process.[45] During femtosecond laser irradiation, free electrons in the metal absorb laser energy, and electron temperature increases rapidly (relative to lattice heat capacity, electron heat capacity is very low).[59,60] Compared with the change in the electron temperature, the change in the lattice temperature is negligible. With the increase of delay time, electrons transfer energy to the lattice through electron-lattice coupling mechanism and then the lattice temperature begins to increase (while the electron temperature begins to decrease), so that the lattice begins to melt, evaporates, and forms plasma. In addition, the kinetic energy of free electrons for the circular laser polarization is higher than that for linear laser polarization.[61,62] After the circularly polarized femtosecond laser irradiates the copper target, the energy of the free electrons in the target will be higher than that obtained under the linearly polarized laser. Therefore, the energy of the electrons transferred to the lattice will be relatively high. The circular polarized laser will produce a higher temperature and density plasma, resulting in a stronger spectral emission.

The mechanism of how the polarization state of laser pulse affects the plasma spectrum has been debated for several years. Previous experimental studies found some differences between plasmas produced by light with different polarization states. For example, Lemos et al. demonstrated experimentally that the expanded velocity of plasma plume produced by femtosecond laser with circular polarization is faster than that with linear polarization.[56] It is known that the femtosecond laser pulses can excite the plasma under the ionization mechanism, which is called “above threshold ionization (ATI)”. The velocity of ionized electrons is determined by the phase of laser electric field at the moment of ionization. When electrons are ionized into free electrons, the energy acquired by electrons is related to the phase of electric field. Corkum et al. pointed out that the electron temperature of an under-dense plasma could be controlled through changing the polarization and wavelength of the ionizing laser pulse.[3] They also calculated the electron energy spectra, respectively, for circularly and linearly polarized laser-induced plasma, discovering that the spectral energy for circularly polarized laser is higher and therefore gives rise to a higher plasma temperature which in turn accelerates the plasma expansion. Arissian et al. investigated the kinetic energy of electrons of plasma created by the circularly and linearly polarized laser pulses with circular and linear polarizations.[40] According to their results, the electron momentum distribution measured after tunnel has been ionized with short pulses, and the associated THz radiation show that the distribution of energy produced by the circularly polarized laser is different from that by the linearly polarized laser. Circularly polarized light leads to higher kinetic energy.[63,64] Therefore, the collisions between electrons and Cu atoms are more frequent in circular polarization condition, creating the denser plasma. To verify the above discussion, we measure the distribution of spectral intensity of Cu plasma, which is excited by laser pulses with different polarizations for different laser energ values (see Fig. 3). The results support them by showing that the density of Cu plasma is higher in the case of circularly polarized laser.

Fig. 3. Variation of Cu (I) lines at (a) 510.55 nm and (b) 521.82 nm with angle of quarter-wave plate for different laser energy.

To further discuss the effect of laser polarization on the spectral emission, we present the variation of Cu (I) line at 521.82 nm with the rotation angle of quarter-wave plate for different laser energy values as shown in Fig. 4. From this figure it follows that the spectral emission in the circular polarization is the highest while that in the linear polarization is the lowest. The Cu atoms will collide with free electrons with relatively high kinetic energy, since electrons can obtain higher energy from the circularly polarized pulses while electrons have lower kinetic energy from the linearly polarized one.

Fig. 4. Evolutions of Cu (I) line at 521.82 nm with the rotated angle of quarter-wave plate for different laser energy values. The laser energies are 0.5 mJ (a), 1.0 mJ (b), 1.5 mJ (c), and 1.9 mJ (d).

The question is how the polarization states affect the kinetic energy of electrons. According to the Mitryukovskiy et al.’s study,[55] for linearly polarized laser pulses, free electrons possess low kinetic energy because they are alternately accelerated and decelerated by laser electric fields in each optical cycle of the pulses. By comparison, electrons are continuously accelerated under the action of circularly polarized pulses, which enables the electrons to collide with more Cu atoms, and increases the emission intensity of breakdown spectrum. To be more specific, a large number of electrons with a kinetic energy of about 2Up are first generated at the end of laser pulse ( is the ponderomotive potential of the electron in laser electric field with linear polarization, with ε0 being the vacuum permittivity, me the electron mass, I the intensity, and ω0 the frequency of the laser field), and the distribution of traverse electron kinetic energy can be estimated through semianalyitical laws (the combination of Newton’s equations for electron motion). At the end of the laser pulse, the energy of most electrons is less than 1 eV in the case of linearly polarized laser. In contrast, the energy is close to 14.6 eV in the the case of circularly polarized laser.

From the above discussion, the electron kinetic energy is influenced by the polarization of light since circularly and linearly polarized laser pulses are subjected to different acceleration and deceleration processes, and the acceleration process can enhance the energy of the charged particles. The fundamental reason of this phenomenon still needs exploring. Bashmakov et al. provided us with guidance in quantum electrodynamical cascades.[65] They investigated the effect of laser polarization on quantum electrodynamical cascading in a standing electromagnetic wave through numerical simulation, indicating that the cascade dynamics is obviously more complicated in the laser field in the case of linearly polarized laser since the electrodynamical field is divided into two qualitatively different parts: the “electric” regions (where the electric field is stronger than the magnetic field) and the “magnetic” regions (where the magnetic field is stronger than the electric field). In the “magnetic” regions, acceleration is an obstacle, while the high-energy electrons are cooled down through emitting photons. On contrary, acceleration is steady in the “electric” region. The volumes of “magnetic” and “electric” region evolve periodically in the case of linearly polarized laser while the volum remains invariable over time in the case of circularly polarized laser. It is noted that acceleration increases the probability of quantum electrodynamical cascading process significantly by turning around the particles transversely to the field and the cascade rate is increased accordingly. Therefore, the number of cascade particles in the case of circularly polarized laser undergoes an exponentially increase while the number of particles in the case of linearly polarized laser grow in time stair-step-like on a logarithmic scale. The energy spectra of electrons under linearly polarized laser pulses show a moderate increase in the particle number accompanied by plasma heating under the action of “electric field”, and it peaks when the whole space is covered by the “electric” region, and become minimal as the “electric” field is replaced completely by the “magnetic” field. Therefore, it can be concluded that the different physics processes of “electric” region and “magnetic” region lead to different acceleration and deceleration processes, which causes electron kinetic energy to be higher (in the case of circularly polarized laser) or lower (in the case of linearly polarized laser) and eventually affects the collisions of electrons with Cu atoms.

As mentioned in the above theory, the energy of free electrons produced by the circularly polarized laser pulse is higher than that by the linearly polarized light, which in turn causes more collisions and thus heats the plasma. We attempt to obtain the plasma temperature and electron density for supporting our propositions. We assum that the laser-induced plasma is in local thermodynamic equilibrium and optically thin, the plasma temperature is calculated by the Boltzmann plot method through using the integrated intensity of the emission line emission of copper plasma: ln(λI / gkAki) = –Ek / kBTe + C,[66,67] where λ is the wavelength, I the line intensity, k the upper level, i the lower level, gk the degeneracy of the upper level, Aki the transition probability, Ek the energy, kB the Boltzmann constant, Te the plasma temperature, and C the constant. In order to obtain the plasma temperature, the Cu (I) 510.55-nm, 515.32-nm, and 521.82-nm lines are selected and used to calculate the plasma temperature. This is because the three spectral lines have all a wide upper level energy range, which can highly improve the accuracy of the experimental result.[68] Relevant parameters are obtained from previously published papers.[45,69] The numerical value of ln(λI/gA) is calculated to estimate the plasma temperature. The obtained data and Ek are plotted, and then the slope of function -(1/kBTe) is obtained by fitting the coordinates of linear points. The typical linear fitting diagram is presented in Fig. 5. Finally, the plasma temperature is obtained by the slope value.

Fig. 5. Typical Boltzmann plots for selected angles with laser energy being 1.9 mJ. The angles are (a) 0°, (b) 45°, (c) 90°, and (d) 135°.

Now consider the electron density. Here, the electron density is calculated by the FWHM of the emission line and Stark broadening: Δλ1/2 = 2ω (Ne / 1016), where Δλ1/2 is the FWHM of emission line, ω the electron collision coefficient, and Ne the electron density. In addition, the line broadening includes natural broadening (ΔλN), Doppler broadening (ΔλD), instrument broadening (Δλinst), Δλ1/2 = Δλobs –ΔλN – ΔλD – Δλinst, with Δλobs being the observed line width in the experiment. In laser-induced plasma, the natural broadening and Doppler broadening are very small compared with the Stark broadening and can be ignored. The line width Δλ1/2 = Δλobs – Δλinst. The observed line width is about 0.04 nm, which is determined by measuring the width of Hg lines emitted from a mercury lamp under low pressure. We use Cu (I) 521.82-nm line to calculate the electron density. The typical Lorentz fitting curves for the observed line width (Δλobs) at different angles of 1/4 wave plate (0, 45, 90, and 135) is shown in Fig. 6.

Fig. 6. Typical Lorentz fitting for selected angles with laser energy being 1.9 mJ and angle being (a) 0°, (b) 45°, (c) 90°, and (d) 135°.

Figure 7 shows the evolution of the plasma temperature and electron density from copper with the angle of quarter-wave plate. The change in the plasma temperature is similar to that in the electron density. As shown in the figure, under circularly polarized femtosecond laser, the temperature and density reach maximum, and they decrease to the minimum under linearly polarized femtosecond laser, while the value in the case of elliptical polarization is between the values of linearly and circularly polarized pulses. This phenomenon supports the previous conclusion that those collisions of electrons are powered by circularly polarized laser pulses and Cu atoms give rise to higher electron density and plasma temperature. These electrons have higher kinetic energy than those produced by linearly polarized laser, which enables them to collide with more Cu atoms, and thus producing more electrons and rasing the plasma temperature. Finally, figure 8 presents the evolution of signal-to-background for Cu (I) 521.82 nm with the angle of quarter-wave plate, with laser enegy being 1.9 mJ. The signal-to-background ratio in the case of circularly polarized laser is higher than that in the case of linearly polarized laser. It indicates that the increase in spectral peak intensity is more than that in background continuous emission. The circularly polarized laser is suitable for element analysis in femtosecond LIBS.

Fig. 7. Evolution of (a) plasma temperature and (b) electron density with angle of quarter-wave plate, with laser energy being 1.9 mJ.
Fig. 8. Evolution of signal-to-background with angle of quarter-wave plate for Cu (I) 521.82 nm, with laser energy being 1.9 mJ.
4. Conclusions

In this work, we demonstrate that the spectral line intensity of femtosecond laser-induced Cu breakdown spectrum strongly depends on laser polarization. The circularly polarized femtosecond laser leads to higher spectral intensity than the linearly polarized femtosecond laser due to the presence of electrons with high kinetic energy created by the circularly polarized laser pulse. The evolutions of two regions (the “electric” one and the “magnetic” one) are different, which gives rise to different acceleration and deceleration mechanism which can account for the difference in electron motion. Those electrons with higher kinetic energy are more likely to collide with Cu atoms and produce higher temperature plasma, causing higher emission intensity. We believe that these findings are conducive to understanding the mechanism of femtosecond LIBS and improving the sensitivity of this technique.

Reference
[1] Miziolek A W Palleschi V Schechter I 2006 Laser-induced breakdown spectroscopy Cambridge Cambridge University Press
[2] Cremers D A Radziemski L J 2006 Handbook of laser-induced breakdown spectroscopy John Wiley & Sons
[3] Corkum P B Burnett N H Brunel F 1989 Phys. Rev. Lett. 62 1259
[4] Liu J Tao H Y Gao X Hao Z Q Lin J Q 2013 Chin. Phys. 22 044206
[5] Yang H X Fu H B Wang H D Jia J W Sigrist M W Dong F Z 2016 Chin. Phys. 25 065201
[6] Yang D P Li S Y Jiang Y F Chen A M Jin M X 2017 Acta Phys. Sin. 66 115201 in Chinese
[7] Yang X Li S Jiang Y Chen A Jin M 2019 Acta Phys. Sin. 68 065201 in Chinese
[8] Wu Z Zhang N Zhu X An L Wang G Tan M 2018 Chin. Phys. 27 077901
[9] Zheng P C Li X J Wang J M Zheng S Zhao H D 2019 Acta Phys. Sin. 68 125202 in Chinese
[10] Wang Y Chen A Li S Sui L Liu D Tian D Jiang Y Jin M 2016 J. Anal. At. Spectrom. 31 497
[11] Gautier C Fichet P Menut D Lacour J L L’Hermite D Dubessy J 2005 Spectrochim. Acta B: At. Spectrosc. 60 265
[12] Sattmann R Sturm V Noll R 1995 J. Phys. D: Appl. Phys. 28 2181
[13] Stratis D N Eland K L Angel S M 2001 Appl. Spectrosc. 55 1297
[14] Benedetti P Cristoforetti G Legnaioli S Palleschi V Pardini L Salvetti A Tognoni E 2005 Spectrochim. Acta B: At. Spectrosc. 60 1392
[15] Zhang D Chen A Wang X Li S Wang Y Sui L Jiang Y Jin M 2017 Opt. Laser Technol. 96 117
[16] Nassef O A Elsayed-Ali H E 2005 Spectrochim. Acta B: At. Spectrosc. 60 1564
[17] Rai V N Singh J P Yueh F Y Cook R L 2003 Laser Part. Beams 21 65
[18] Harilal S Tillack M O’shay B Bindhu C Najmabadi F 2004 Phys. Rev. 69 026413
[19] Pandey P K Thareja R K 2013 Phys. Plasmas 20 022117
[20] Singh K S Sharma A K 2016 Phys. Plasmas 23 122104
[21] Popov A M Colao F Fantoni R 2010 J. Anal. At. Spectrom. 25 837
[22] Popov A M Colao F Fantoni R 2009 J. Anal. At. Spectrom. 24 602
[23] Shen X K Sun J Ling H Lu Y F 2007 Appl. Phys. Lett. 91 081501
[24] Guo L Hao Z Shen M Xiong W He X Xie Z Gao M Li X Zeng X Lu Y 2013 Opt. Express 21 18188
[25] Wang Q Chen A Zhang D Wang Y Sui L Li S Jiang Y Jin M 2018 Phys. Plasmas 25 073301
[26] De Giacomo A Gaudiuso R Koral C Dell’Aglio M De Pascale O 2013 Anal. Chem. 85 10180
[27] Dell’Aglio M Alrifai R De Giacomo A 2018 Spectrochim. Acta Part. B: At. Spectrosc. 148 105
[28] Freeman J R Harilal S S Diwakar P K Verhoff B Hassanein A 2013 Spectrochim. Acta B: At. Spectrosc. 87 43
[29] Ahamer C M Pedarnig J D 2018 Spectrochim. Acta B: At. Spectrosc. 148 23
[30] Li K Zhou W Shen Q Ren Z Peng B 2010 J. Anal. At. Spectrom. 25 1475
[31] Hou Z Wang Z Liu J Ni W Li Z 2014 Opt. Express 22 12909
[32] Wang Y Jiang Y He X Chen Y Li R 2018 Spectrochim. Acta B: At. Spectrosc. 150 9
[33] He X Chen B Chen Y Li R Wang F 2018 J. Anal. At. Spectrom. 33 2203
[34] Zhang D Chen A Wang Q Wang Y Qi H Li S Jiang Y Jin M 2018 Phys. Plasmas 25 083305
[35] Guo J Wang T Shao J Chen A Jin M 2018 J. Anal. At. Spectrom. 33 2116
[36] Sanginés R Sobral H Alvarez-Zauco E 2012 Spectrochim. Acta B: At. Spectrosc. 68 40
[37] Eschlböck-Fuchs S Haslinger M J Hinterreiter A Kolmhofer P Huber N Rössler R Heitz J Pedarnig J D 2013 Spectrochim. Acta B: At. Spectrosc. 87 36
[38] Wang Y Chen A Wang Q Zhang D Sui L Li S Jiang Y Jin M 2019 Plasma Sci. Technol. 21 034013
[39] Wang Y Chen A Sui L Li S Wang X Jiang Y Huang X Jin M 2016 J. Anal. At. Spectrom. 31 1974
[40] Arissian L Mirell D Yeak J Rostami S Diels J C 2012 Multiphoton Processes and Attosecond Physics Berlin, Heidelberg
[41] Lee Y Mao X Chan G C Y Gonzalez J Russo R E Zorba V 2018 J. Anal. At. Spectrom. 33 1875
[42] Harilal S S Yeak J Brumfield B E Suter J D Phillips M C 2016 J. Anal. At. Spectrom. 31 1192
[43] Hou H Cheng L Richardson T Chen G Doeff M Zheng R Russo R Zorba V 2015 J. Anal. At. Spectrom. 30 2295
[44] Hou H Mao X Zorba V Russo R E 2017 Anal. Chem. 89 7750
[45] Wang Q Chen A Qi H Li S Jiang Y Jin M 2020 Opt. & Laser Technol. 121 105773
[46] Nakimana A Tao H Y Hao Z Q Sun C K Gao X Lin J Q 2013 Chin. Phys. 22 014209
[47] Semerok A Dutouquet C 2004 Thin Solid Films 453 501
[48] Angel S M Stratis D N Eland K L Lai T Berg M A Gold D M 2001 Fresenius�?Journal Analytical Chemistry 369 320
[49] Yang F Jiang L Wang S Cao Z Liu L Wang M Lu Y 2017 Opt. Laser Technol. 93 194
[50] Pinon V Fotakis C Nicolas G Anglos D 2008 Spectrochim. Acta B: At. Spectrosc. 63 1006
[51] Pinon V Anglos D 2009 Spectrochim. Acta B: At. Spectrosc. 64 950
[52] Cao Z Jiang L Wang S Wang M Liu L Yang F Lu Y 2018 Spectrochim. Acta B: At. Spectrosc. 141 63
[53] Kim J Na S 2007 Opt. Laser Technol. 39 1443
[54] Wang Y Chen A Wang Q Zhang D Li S Jiang Y Jin M 2020 Opt. Laser Technol. 122 105887
[55] Mitryukovskiy S Liu Y Ding P Houard A Couairon A Mysyrowicz A 2015 Phys. Rev. Lett. 114 063003
[56] Lemos N Grismayer T Cardoso L Geada J Figueira G Dias J M 2013 Phys. Plasmas 20 103109
[57] Wang Q Chen A Xu W Li S Jiang Y Jin M 2019 J. Anal. At. Spectrom. 34 1242
[58] Qi H Li S Qi Y Chen A Hu Z Huang X Jin M Ding D 2014 J. Anal. At. Spectrom. 29 1105
[59] Wang T Guo J Shao J Wang D Chen A Jin M 2015 Phys. Plasmas 22 033106
[60] Li S Li S Zhang F Tian D Li H Liu D Jiang Y Chen A Jin M 2015 Appl. Surf. Sci. 355 681
[61] Shi Y Chen A Jiang Y Li S Jin M 2016 Opt. Commun. 367 174
[62] Li S Y Li S C Sui L Z Jiang Y F Chen A M Jin M X 2016 Phys. Rev. 93 013405
[63] Smeenk C T Arissian L Zhou B Mysyrowicz A Villeneuve D M Staudte A Corkum P B 2011 Phys. Rev. Lett. 106 193002
[64] Zhou B Houard A Liu Y Prade B Mysyrowicz A Couairon A Mora P Smeenk C Arissian L Corkum P 2011 Phys. Rev. Lett. 106 255002
[65] Bashmakov V F Nerush E N Kostyukov I Y Fedotov A M Narozhny N B 2014 Phys. Plasmas 21 013105
[66] Wang Y Chen A Wang Q Sui L Ke D Cao S Li S Jiang Y Jin M 2018 Phys. Plasmas 25 033302
[67] Wang Q Chen A Wang Y Sui L Li S Jin M 2018 J. Anal. At. Spectrom. 33 1154
[68] Zorba V Mao X Russo R E 2015 Spectrochim. Acta B: At. Spectrosc. 113 37
[69] Wang Y Chen A Zhang D Wang Q Li S Jiang Y Jin M 2020 Phys. Plasmas 27 023507